Far-field interaction of focused relativistic electron beams in electron energy loss 

spectroscopy of nanoscopic platelets 
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A quantum mechanical scattering theory for relativistic, highly focused electron beams near 
nanoscopic platelets is presented, revealing a new excitation mechanism due to the electron wave 
scattering from the platelet edges. Radiative electromagnetic excitations within the light cone are 
shown to arise, allowed by the breakdown of momentum conservation along the beam axis in the 
inelastic scattering process. Calculated for metallic (silver and gold) and insulating (Si02 and MgO) 
nanoplatelets, new radiative features are revealed above the main surface plasmon-polariton peak, 
and dramatic enhancements in the electron energy loss probability at gaps of the 'classical' spectra, 
are found. The corresponding radiation should be detectable in the vacuum far-field zone, with 
e-beams exploited as sensitive 'tip-detectors' of electronically excited nanostructures. 

PACS number(s): 79.20.Uv, 78.67.Bf, 73.20.Mf, 41.60.-m ' 
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I. INTRODUCTION 

A powerful technique for investigating electromag- 
netic (EM) field distribution around nanostructures is 
provided by very fast (relativistic) electron beams (e- 
beams), with typical lateral resolution on an atomic 
scale, available in scanning transmission electron micro- 
scopes (STEM)ii£iM£i&. As discussed previously, when 
the e-beam is restricted to the vacuum near a selected 
nanoparticle-' 2 , its EM interaction with surface plas- 
mons or surface plasmon-polaritons (SPPs)£ is reminis- 
cent of the near-field interaction^ of subwavelength opti- 
cal probes. Several works have recently studied realiza- 
tions of Cherenkov radiation excitation within various 
dielectric media by e-beams moving in near-field vac- 
uum zone a 8 ' 9 ' 10 ! 11 . In all the latter works the energy 
loss processes were described within a simple classical 
model in which the fast electron was assumed to move 
with a constant velocity along a straight line trajectory 
near a finite dielectric medium such that the energy loss 
intensity could be obtained from the force exerted on 
the electron due to its self-induced electric field through 
the nearby dielectric medium. The great simplification 
achieved by this approach amounts to reducing the full 
scattering problem at hand to a problem of finding the 
EM field induced by the e-beam in the vacuum around 
the dielectric medium. The resulting EM field could 
include Cherenkov-like radiative components around the 
e-beam which were restricted, however, to propagation 
within the interior of the dielectric medium. 

In this paper we present a quantum mechanical theory 
for the inelastic scattering of a relativistic highly focused 
e-beam traveling near nanoparticles in a 'non-touching' 
aloof configuration^ 2 -. We show that the electron wave 
scattering by nanoparticle edges along the beam axis 
switches on Chcrcnkov-likc radiation channels which ex- 
tend into the vacuum away from the nanoparticle. The 
resulting far-field coupling between the electron and the 
nanoparticle is found to dramatically enhance various ra- 



FIG. 1: Schematic illustration of the e-beam configuration 
with respect to the rectangular platelet used in our calcu- 
lations. Another configuration, used in Ref.— , where the 
e-beam is parallel to an edge of the platelet is also shown. 



diativc channels in the loss spectrum. 

To illustrate our main points we consider here a simple 
model (see Fig.l) where the e-beam is propagated in the 
vacuum along a wide face of a rectangular nanoplatelet 
(oriented, e.g., in the x — y plane), and a surface or guided 
wave induced by the electron is propagated with a wave 
number k x along the beam axis. The spatially sensi- 
tive nature of the corresponding electron energy loss pro- 
cess arises from the exponential dependence, e~ 2K b , of 
the EM interaction between the e-beam and the platelet 
on the impact parameter b. The extinction coefficient, 

K* = \J K 2 - (u>/c) 2 , with K 2 = k 2 + fc 2 , determines 
the tail of the evanescent field in the vacuum for values 
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of k outside the light-cone, i.e. for K > uj/c. Inside the 
light-cone, i.e. for K < lo/c , K* is purely imaginary and 
the corresponding interaction becomes spatially oscillat- 
ing, allowing the electron to exchange photons with the 
particle far away into the vacuum. This striking mech- 
anism has been overlooked in the recent literature of 
STEM-electron energy loss spectroscopy (EELS), since 
the excitation by an electron moving in the vacuum with 
a classical velocity v, has been restricted to a constant 
longitudinal wavenumber k x = tu/v > lo/c, implying EM 
coupling to the nanoparticle which is restricted to the 
evanescent tail near the surface. 

Our model calculation is applied to two types of 
nanoscopic platelets, conducting platelets made of sil- 
ver or gold, and dielectric platelets made of insulators 
such as silica or magnesia. For both types of nanopar- 
ticles we find significant loss signals in the low energy 
range of the spectrum, where the electron-hole excitation 
probability is either zero (for the insulator) or very small 
(for the metals), exhibiting far-field (radiative) charac- 
teristics. In particular, specific SPP modes of the silver 
platelet, which penetrate into the light cone, can be ex- 
cited by the external e-beam, leading to new features in 
the EEL spectrum which decay weakly with the beam- 
platelet distance. 



position vector, {s} is a collective symbol for the position 



vectors of the platelet charges, and m — rriQ / y 1 — (v /c) 
the dynamic electron mass. 

To first order of the perturbation theory with respect 
to the EM interaction Hamiltonian, Hem , the probabil- 
ity for the e-beam to go, during the time interval r, from 
initial to final eigen states when the platelet initial state 
is the ground state is: 
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II. MODEL AND FORMULATION 

Following Ref4, the focused e-beam is described here 
as a one-dimensional wave, propagating along the x-axis, 
while in the transverse (y — z) directions it is described 
by a wave-packet localized within a smoothly converg- 
ing cross section along the beam axis, whose shape is as- 
sumed to be squared for the sake of simplicity. The corre- 
sponding Green's function for the noninteracting focused 
e-beam may be therefore written in the general form: 



2L 

Xt tr (V,z]x) x^ tr {y\z';x') 



(1) 



where p x = hq x is the longitudinal (along the beam axis) 
electron momentum, £ Px , Ptr = \/p x c 2 + rri^c A + c 2 p\ r its 
total relativistic energy eigenvalue, with mo the electron 
rest mass, ~p tr = ftqtr ( ~Qtr = (Qy,Qz) ) its transverse 
momentum, and e. lqmX X~q tl . {V-i z 'i x ) t ne corresponding e- 
beam eigenfunction (see Appendix A). 

The electromagnetic (EM) interaction between the e- 
beam and the platelet may be described effectively by 
the Hamiltonian: 



Hem ( r>, {*}) « -e* 0*. W) " e— A x (V, {s}) (2) 

mc 

where $ and A x are the scalar and x-component of the 
EM four-vector potential respectively, r is the electron 



where the the sum is over the platelet final states af 
, and : A(e,e';r) = -MM^Mzl 
oo , 



In the limit 



when: t — » oo , |A(e,e';r)| — ► 2tt hr 5 [e ~ e') — 
2t Re J °° dt exp [it (e — e') /ft] , and so the rate of change 
of scattering probability of the e-beam R e u_>t-\ = 



K 



-V 

is given by: 
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where the inclusion of all terms under the real part sym- 
bol is justified by the reality of the total expression writ- 
ten in the last five rows within the curly brackets. 

Using the relations: 
e<«*-f x H p J M (?), 
(— e) $ (V) — (V) , the rate of change of 

probability for the scattering of the e-beam can be 
rewritten in the form: 



with: H E ° M ("r ) 
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rearranging the integrations we find that: 
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where 



ton Green's function D v ^ ( 
(«-» ct, a;, y, z), as: 



em j J/' j t) H P i M {x, y, z; 0)^> 
D ,o (?',?;t) + £-D QA (1», 
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Aq x = (g< - qf) « («/«) + ft (g/ r ) 2 - (<z t V) 2 /2nw 

(6) 

is the longitudinal momentum transfer of the e-beam 
(see Appendix A), ~q\ r and ~~q{ r the e-beam asymptotic 
transverse momenta, initial and final respectively, and 

Huj = \e v i -ffi — st -*i ) its energy loss. Note that the 

\ fx' 1 tr Pxi Q tr J 

width /3 _1 of the Gaussian distribution function, is intro- 
duced in Eq. |(5J) to account for the high transverse-energy 
cutoff caused to the e-beam by the objective aperture. It 
is related to the length L of the region around the beam 
focal plane used in our model as a normalization factor 
for the electron wave functions. 

The interaction potential, Hem (x, y, z) , between the 
platelet and an external electron at (x, y, z) is nearly in- 
dependent of x for |x| <C a*, and decays to zero at least 
as quickly as 1/x 2 for \x\ > a* (see, e.g., RefJ^). Under 
these circumstances the limits of the integrations over x 
and x' in the above expression may be set at —a* and 
a*, rather than at — L and L. The correlation function 

Hjsm y'l 2 '' ^™ ( x > y> 2 i 0)\ can be expressed in 
terms of the relevant components of the 4-tensor pho- 



;*), v,[i = 0,1,2,3 



For the sake of simplicity we may assume trans- 
lational invariance of the platelet dielectric properties 
in the x — y plane, that is take: D V!fl , r ;t) — 
Du.fi (x — x', y — y', z' , z; t). For an impact parameter b 
smaller than the platelet sides along the x and y axes 
(i.e. b <§; 2a* , 2b* ) this assumption may be justified, 
though it is inconsistent with the breakdown of momen- 
tum conservation in the beam-platelet scattering event 
considered here (see a more detailed discussion below). 

Substituting into the above expression for R e n_^f\ and 
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Jdzfdye lk " v x*-,f {y,z;x)xtt (y, z ;x) 

1 tr 



Dd,Pl (k x ,k y ,uj;z',z)} 

where D v ^ iV ^ (k x , k y ,u>; z', z) 

= f^dte^DP^P* (k x ,k y ;z',z;t), 

D^(k x ,k y ;z',z;t) = (7) 

A),o (k x , k y ; z', z; t) + %^D ,i {k x , k v ;z', z; t) + 
^Di,o (k x ,k y ; z', z; t) + ^,J)i,i (k x , k y ; z', z; t) 

and D ll ^{k x ,k y ; z 1 , z;t) is the spatial Fourier trans- 
form of D v ll (x — x' , y — y' , z' , z; t) with wavevector K — 
[k x , k y ). 

Now, the 4-tensor photon propagator in the vacuum 
(i.e. at z, z' < ) has the forrn^: 

Di/ t n {k x , k y , lu; z , z) 



Vu 



2nK* 



's v ,^e- K *\''-'\ - r Vj/i (k x , k y , u) e K ^ z ' +z ) 



in which the relevant part is associated only with the 
second term within the square brackets (i.e. that associ- 
ated with the image potential of the e-beam) . Using this 
expression and recalling that 1/ \K*\ is typically much 
larger than the beam transverse dimension, so that the 
extreme confinement of the e-beam wave functions X^-s 

under the integrals over z and z' restrict their values to 
near z' = z = —b ,we have: 



a narrow region 



R, 
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IjS) I 

J* ( 1 tr) Qtr'i kyi K*; a;' J 

2L J-a* " XC 
J ( Q tr ) 9tr i ' ; X 



where: 



r f ' 1 (k x , k y ,uj) = r , (k x ,k y ,uj) + ^-r ,i (k x ,k y ,cu) + 



rri.o (k Xl k y ,uj) + ( J' c) V i,i (k x ,k y ,Lu) 



(8) 
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and: 

j(y{ r ,qi;k y ,K*;x) (9) 
dz J dye- lk y v X *-ff r (y, z; x) Xqlr {y, z\ x) 
Finally, denoting: 

l{lH-T tr -k y ,K*-{Aq x -k x )) (10) 

= lL dxeliAq *~ k * )X J (Jill Ql; k v ,K*; 
the scattering rate is rewritten as: 



constant, that is equivalent to a S- function in the corre- 
sponding real-space transverse coordinates, (2) the con- 
tribution of the transverse energy to the longitudinal mo- 
mentum transfer Aq x (see Eq.©) can be neglected, and 
(3) the effective particle size, a* , appearing as an integra- 
tion limit along the beam axis, is infinite. Assumption 
(3), in conjunction with (1), yields the conservation of 
longitudinal momentum, i.e. Aq x — k x = 0, which to- 
gether with assumption (2) imposes the fixed condition 
k x = (uj/v). 

It is interesting to note that usually assumption (2) 
is not strictly satisfied since the contribution of the 



transverse energy to Aq x : h 



(it 



/2mv 



R, 



e(i- 



— — / dk T / dk 



Im 



2e 2 
(k x , ky,uj) 



(11) 



-2K*b 



K* 



E 



~Qtr\ 



n 1 -k K*- 



(Ag* 



qtrAqtr/ (mv/h) ~ ±q 2 r /q l x , can be as large in magni- 
tude as (w/v). As an example, at ftw ~ 10 eV, (co/v) 
~ 0.05 run -1 , whereas the transverse beam-wavenumber 
uncertainty, |Ag tr | ~ q tr ~ 2tt/Z (with a typical value of 
I ~ 0.6 nm for the beam radius) is 10 nm -1 , so that for 

= 100 keV , where q x ~ 1500 nm" 1 , q 2 tr jq x 
l 



£ 

nm 



0.07 



III. THE 'CLASSICAL' APPROXIMATION AND 
BEYOND 

The theory developed in the previous section can be 
further simplified without losing its main physical con- 
tent by employing several approximations. In the long 
wavelengths limit discussed in Ref* 1 ^ we find that (see 
Appendix B): 



Im 



Im • 



-2K*b„ 



/K* 



{ [(K */k 2 ) f e + ([v/cf - (co/ck) 2 



fo/K* 
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where 

fe= ' 

At = 

Dt = 



e 2 K* 2 



Q 2 ) /D+D-, f = (K* 2 - Q 2 ) /D+D a 



eK* + Q tanh (Qc*) , D~ = eK 
K* + Q tanh (Qc*) , D~ = K* -\ 



+ Q coth (Qc*) 
Qcoth (Qc*) 



Q = y 'K 2 — (uj j c) 2 e (uS) , and e (to) is the local bulk 
dielectric function of the platelet. In the limit of a semi- 
infinite medium the resulting expression reduces (see Ap- 
pendix B) to the surface dielectric response function 
obtained in Ref, 14 by using Maxwell's equations with 
macroscopic boundary conditions. 

The standard classical approximation for the loss 
function 14 is obtained from Eq. lfTTj) by making the fol- 
lowing assumptions: (1) the e-beam transverse momen- 
tum distribution function J 



In the present paper we focus on the most interesting 
violation of the 'classical' approximation outlined above, 
allowing a* to be a finite length, which reflects an effec- 
tive range of the actual beam-particle interaction along 
the beam axis. Consequently the longitudinal momen- 
tum distribution around Aq x — k x — 0, defined by the 
integral in Eq. (jl0p . is smeared and many wavenumbers 
k x inside the light-cone start contributing to the loss rate, 
Eq.GD. 

The condition for the smearing to be significant is 
Tt/a* > (uj/v), so that typically for frequencies u> in 
the visible range, a* should be smaller than 200 nm. 
Nanoplatelets of those lengths should dramatically en- 
hance radiative excitations by the e-beam, previously 
overlooked in the literature; see e.g. Ref.— where it was 
(12^rgued that recoil effects in STEM should be negligible 
^ for valence electron excitations. Recoil is only a clas- 
sical remnant of the present effect and of less general 
appearance. In particular, it vanishes for large media, 
such as the porous film investigated in Ref. 8 , for which 
(if made sufficiently thin) the quantum mechanical mo- 
mentum uncertainty along the e-beam axis remains sig- 
nificant. 

It should be stressed that, for the sake of simplic- 
ity, the platelet dielectric response is calculated by as- 
suming its wide faces to be infinite. A fully consis- 
tent treatment of the breakdown of translation invari- 
ance is expected, however, to further enhance all radia- 
tive channels. The calculation of the kinematical factor 
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111 



y tri 



k K*- 



q \ T \ k y , K*; (Aq x — k x )j , responsible for the lon- 
gitudinal momentum uncertainty in our model, from the 
integral in Eq. lTlT))) . could generate artificial oscillations 
by the sharp cutoff of the integral at x = ±a*. To avoid 
such oscillations we use an equivalent Gaussian distribu- 
tion function in our actual calculations. The correspond- 
ing smooth cutoff is, in fact, more realistic than that ap- 
pearing in Eq. (|10|) since it arises from the attenuation of 
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the e-beam-platelet interaction at \x\ values larger than 
a*. In any event, the exact form of the corresponding 
distribution function is of no great importance for the 
main purpose of our present paper. 



IV. RESULTS AND DISCUSSIONS 



A. Silver and Gold Nanoplatelets 



As a first example we calculate the EEL function of 
a 100 nm long silver and gold platelets for an external 
100 keV e-beam at various impact parameters (see Figs. 2 
and 3). To analyze the various SPP resonances one may 
consider the zeros of the denominator of the extraordi- 
nary wave amplitude f e in Eq. (|12p in the complex K- 
plane. With the experimental optical dielectric function, 
e{uS), for silver— the resulting dispersion relation (inset, 
Fig. (2)) exhibits a rather flat branch of ui(ReK) inside 
the light-cone, which can be attributed to radiative SPP, 
seen as a mirror image of the usual non-radiative SP dis- 
persion curve with respect to the light-line. The sector of 
w(Re K) connecting the two branches across the light line 
has a vanishing negative slope, where IraK(cu) oc Ime(u) 
has a sharp peak. The sharp dip in the EEL spectrum 
just above the classical SP frequency (at 3.8 eV) reflects 
these closely related features. 

At slightly higher frequencies the EEL signals exhibit 
a pronounced rise due to the enhanced SPP density of 
states associated with the flat radiative SPP branch. 
These peculiar features are missing in the loss spectra 
of the gold platelet, shown in Fig. 3. 

The EEL intensity in this spectral region exhibits at- 
tenuation with increasing impact parameter significantly 
weaker than the corresponding attenuation of the main 
SP peak calculated in the classical limit. The radiative 
nature of the beam-particle coupling shown in Figs. (2,3) 
is even more pronounced in the low energy region below 
the main SP peak, where the classically calculated sig- 
nal drops to very small values. Here our calculated EEL 
function exhibits a pronounced broad band with linearly 
increasing intensity for increasing frequency and almost 
no attenuation with increasing impact parameter. These 
features are due to the fact that the loss signal well be- 
low the main SP frequency is dominated by the contri- 
bution from the ordinary wave amplitude f , appearing 
in Eq.(fT2|). which is singularly enhanced near the light 
line (where K* —* ), and thus reflecting the nearly 
pure (transverse) photonic nature of the excitations by 
the e-beam in this 'classically forbidden' region. 

The results of our calculations may be compared to 
the experimental data reported in Refpii for silver and 
gold nano rods and ellipsoids. Fig. 4 shows our calculated 
EEL spectra for three silver platelets with c* — 15 nm 
and a* — 10, 15,30 nm at impact parameter b = 10 nm. 
The shown curves may be compared to the spectrum in 
RefJi obtained for a silver ellipsoid with a long half-axis 



FIG. 2: EEL spectra (solid lines) of a 100 keV e-beam propa- 
gating parallel to the a;-axis of a rectangular Ag platelet (with 
half sides: a* — 50 nm along x , and c* = 10 nm along z , see 
Fig.l) at impact parameters b = 10, 20,40 nm above its wide 
x — y face. The experimental optical dielectric function, e(ui), 
for silver— has been exploited. Dashed lines represent spec- 
tra calculated by the classical theory. Inset: surface plsmon 
polariton (SPP) dispersion curves, u) (ReK) , u> (Im.fr) in the 
complex K -plane for silver. The indicated values of K and u 
are normalized by K n = u„/c, and w n = 10 eV, respectively. 



FIG. 3: The same as Fig. 2 for a platelet made of gold (Au). 
Note the absence of the sharp dips appearing just above the 
main plamon peaks in the corresponding Ag spectra (Fig. 2). 
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FIG. 4: The same as Fig. 2 for three Ag platelets with half 
sides along the beam axis a* = 10, 15, 30 nm, and half width 
c* = 15 nm, and with the e-beam at an impact parameter 
b = 10 nm. The dashed line represents the corresponding 
classical spectrum (i.e. for a* — > oo). 



(~ 30 nm ) and two short half-axes (~ 15 nm) at impact 
parameter ~ 10 nm above the ellipsoid wide face. 

The two lower curves, particularly those corresponding 
to a* = 10 nm, exhibit good agreement with the relevant 
experimental data. Specifically, in addition to the very 
good agreement of the calculated main plasmon peak po- 
sition (~ 3.45 eV) with the experimental one, the intensi- 
ties ratio (~ 2) between the main plasmon peak and the 
high energy broad peak, and the extent and magnitude 
of the low energy tail shown in Fig. 4, are seen to agree 
pretty well with the corresponding experimental results. 
In contrast a large intensities ratio (~ 8) and a very small 
low energy tail characterize the classical curve shown in 
Fig. 4, both indicating the importance of the quantum 
effects predicted by our theory. Note that an important 
feature of our calculated spectra, the large dip just above 
the main plasmon peak, which is missing in the experi- 
mental data, is shown to develop only at relatively large 
values of a* (i.e. for a* > 20 nm ). 



B. Insulating Nanoplatelets 

The situation in the forbidden energy gap region of 
semiconductors and insulators is in a sense an extreme 
case of the effect demonstrated in the low energy region 
of Fig. 2: The EEL spectra shown in Fig. 5 are calculated 
for an external 100 keV e-beam, propagating parallel to 
the x — y face of a 100 nm long SiC>2 platelet with half 
thickness c* = 50 nm, at different impact parameters b. 



FIG. 5: EEL spectra (solid lines) of a 100 keV e-beam prop- 
agating parallel to the x-axis of a rectangular Si02 platelet 
at distances b — 2 and 8 nm above its wide (x — y) face. 
The platelet half-sides along the x, and z axes are: a* — 50 
nm, and c* = 50 nm respectively. The corresponding spec- 
tra (dashed lines) obtained from the classical theory are also 
shown for comparison. Note the close similarity of the clas- 
sical spectrum for b = 8 nm with the one obtained in RefJ^ 
for the same e-beam velocity at nearly the same impact pa- 
rameter parallel to a sharp wedge (see Fig.l). 



The spectra reveal a pronounced double-peak structure 
within the forbidden gap region, which does not decay 
with increasing b values. Strictly speaking, this struc- 
ture reduces to a single broad peak for platelets of widths 
c* < 10 nm, reflecting a finite-size effect. Similarly to the 
situation with the silver and gold platelets well below the 
main SP peak, the strong radiative nature of this feature 
arises from the ordinary wave amplitude f a , correspond- 
ing to the excitation of purely transverse EM waves, po- 
larized within the x — y plane, which totally dominates 
the loss signal in the forbidden gap region. 

The spectra shown in Fig. 5 may be compared to the 
results reported in Ref.— for an electron moving paral- 
lel to a 90° SiC>2 wedge at a distance of 8.5 nm (sec 
Fig.(l)). The pronounced radiative broad band within 
the gap region, obtained in our calculation, dramatically 
contrasts the vanishing loss signal shown there in Fig. (4) 
for an electron beam with the same velocity ( v = 0.54c) 
and nearly the same impact parameter. The lack of far- 
field coupling in the latter theoretical approach restricted 
the fast external e-beam to excitation of EM waves con- 
fined within the dielectric medium 8 , similar to ordinary 
waveguide modes which can develop within a thin SiC>2 
slab in the forbidden gap region where Re e (u>) w 2 , 
and Ime (u>) — > 0. For an ideal planar geometry (as as- 
sumed in our calculation of the dielectric response func- 
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tion r yK,u>J), the corresponding waveguide modes ap- 
pear as extremely narrow resonances which can not be ex- 
cited by an e-beam with Aq x values outside the light cone 
due to the vanishingly small dielectric damping, Ime (u>). 

Such radiation excitations become possible for the non- 
planar geometries studied in Refs.— ,— even under the 
rigid e-beam trajectory approximation exploited there 
(but only above a threshold beam energy considerably 
higher than 100 keV) due to the translational symmetry- 
broken dielectric media considered in their calculations. 
Yet, the corresponding Cherenkov-like channels remain 
fundamentally different from the ones we propose: The 
opening of scattering channels with wave numbers inside 
the light cone allows coupling of the e-beam to the con- 
tinuum of EM modes which are extended into the vac- 
uum perpendicular to the platelet wide face. The relative 
strength of the present radiative mechanism may be fur- 
ther appreciated by noting the calculated spectra near 
a sharp SiC>2 wedge in Ref where in spite of the geo- 
metrical enhancement of near-field Cherenkov coupling, 
beam energies far above 100 keV were needed there for 
'switching on' such channels. 

Finally, it is instructive to compare our predicted loss 
spectrum of an external 100 keV e-beam propagating 
above a MgO platelet with half-sides a* = c* = 50 nm at 
an impact parameter 6=2 nm (see Fig. 6) to the exper- 
imental data reported in RefJ^ for a MgO smoke cube 
of 100 nm size. The overall agreement is good, includ- 
ing the occurrence, in both the calculated spectrum and 
the experimental data, of a broad, nonvanishing signal 
within the forbidden gap region, which is missing in the 
classically calculated spectrum. In this gap region the 

calculated spectrum exhibits a smooth oscillatory 
structure associated with the multiple reflection of the 
generated radiation between the two parallel faces of the 
platelet perpendicular to the z-axis. This finite-size ef- 
fect is peculiar to the far-field radiative modes found in 
the present paper for platelets confined in the direction 
along the e-beam axis, and is different from (though re- 
lated to) the extremely sharp resonances associated with 
the waveguide modes developed in an 'ideal' (i.e. wide 
laterally) planar dielectric thin film. Thus, the classical 
approach applied to such an 'ideal' film yields usually 
(i.e. except for extremely rare coincidences of the loss 
energy with the resonant frequencies) null loss intensity, 
whereas in our quantum calculations the continuous win- 
dow of wavenumbers inside the light-cone removes the 
stringent resonant conditions and allows the appearance 
of a significant loss intensity in the entire gap region. It 
is interesting to note that the average calculated signal 
inside the gap region increases smoothly with increasing 
frequency from zero up to nearly the intcrband threshold 
where its intensity relative to the loss main peaks (at ~ 
14 and 20 eV) is about 1/6. This ratio is remarkably 
close to the corresponding relative intensity observed ex- 
perimentally in Refii£. 



FIG. 6: EEL spectra (solid lines) of a 100 keV e-beam prop- 
agating parallel to the a;-axis of rectangular MgO platelets at 
a distance b — 2 nm above their wide (x — y) faces. The 
platelets half sides along the x and z axes are a* — 50 nm, 
and c* — 50, 100 nm, respectively. The corresponding classi- 
cal (a* —> oo ) results (dashed lines) are also shown. Note the 
finite size oscillations of the calculated loss signal inside the 
forbidden energy gap with a period roughly proportianal to 
1/c*. 



V. CONCLUSION 

Applying a quantum-mechanical approach to the scat- 
tering problem of highly focused relativistic e-beams 
near nanoplatelets, we have shown that Cherenkov-like 
radiation of STEM e-beams, discussed recently in the 
literatur o 8 ' 9 ' 10 , has a much broader scope than originally 
presented. Dramatic enhancements of radiative chan- 
nels arise from the breakdown of momentum conserva- 
tion along the e-beam axis in the inelastic process due to 
scattering of the electron wave by the nanoparticle edges. 
Further enhancements, realized due to the extreme lat- 
eral confinement of the e-beam and its associated trans- 
verse momentum uncertainty^, have not been considered 
in detail here. The radiation predicted to be emitted from 
both conducting and insulating nanoplatelets can be gen- 
erated at impact parameters larger than the evanescent 
tail of the excited surface EM modes due to the oscilla- 
tory distance dependence of the electron-platelet interac- 
tion for momentum transfers within the light-cone. Con- 
sequently, this radiation should have a significant prop- 
agation component perpendicular to its main direction 
along the e-beam axis. 

Large deviations from the classical EEL signal are 
found to persist also at small impact parameters, which 
can be readily tested experimentally. The results of our 
calculations for silver platelets seem to agree pretty well 
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with the experimental data reported in RefJ^ for silver 
nano ellipsoids. Furthermore, experimental observation 
of loss signals within the forbidden energy gap of MgO 
cubes of 100 nm size by Aizpurua et alJ^ seems as well 
to support our main prediction. 



VI. APPENDIX A 

In this appendix we specialize our general model of 
the focused e-beam to allow a more detailed discussion 
of some aspects of EELS experiments in STEM pertinent 
to the subject under study in this paper. We employ the 
relativistic Schrodinger's wave equation: 



where 



dy 2 



dz 2 



dx 2 



£ ip(x,y,z) , 
e 2 /H 2 c 2 - m 2 c 2 /h 2 



subject to the boundary conditions: 



ip (x, y,z) = , for I (x) > y > (Al) 
and for I (x) — b > z > —b 

with : I (x) = Iq + a\x\ , aCl 

Due to the small converging angle a one may invoke 
the Born-Oppenheimer approximation: ip(x,y,z) = 
if (x) x (y, z; x), in which the crossed derivatives 

0^x(y,z;x) , g^x(y,z;x) are neglected, and the wave 
equation takes the approximate form: 



x(y, z ;x) 



dy 2 



dz 



x(y,z;x) 



1 d 2 

f (x) dx 2 



tp(x) 



subject to the boundary conditions, Eq. (TT3]). So- 
lutions for the "slow" motion wave equation satisfy- 
ing these boundary conditions are: Xn y ,n z (j/> z 'i x ) — 
j^j sin [q y (x) (y - I {x))\ sin [q z (x)(z + b-l (x))] where: 

q y ,z (x) = Tm y<z /l (x) , n VtZ = 1, 2, .... 

The resulting equation for the "fast" motion is: 



ip(x) 
where: 



e <p(x) 



(A2) 



v lt ( x ) = ll (x) + ql (x) 



Thus, to first order in perturbation theory with respect 
to v 2 -* (x), the energy eigenvalues of an electron 'trapped' 
by the EM lenses inside the conic beam region are given 
by: 



h 2 c 2 e 2 



2 4 



(A3) 



2 2, 2 4, 2 fc 2 / /;\2 

p x c +m Q c + c h (irn/l) 



Vs. 



hq x is the e-beam main (longitudinal) mo- 



mentum, and Py^ z = hmiy^/l , with: 



1,2, ...,l = 



ylo ('o + a L), its transverse momentum components in 
the free propagation zone outside the EM focusing do- 
main. 

This model of the e-beam is, of course, a drastic simpli- 
fication of the actual focused beam in STEM. In particu- 
lar the ideally reflecting boundary conditions, Ea. (fT3)) . 
can not be strictly realized under the smoothly vary- 
ing field generated in space by the focusing EM lenses. 
The results of our analysis here are not expected to be 
very sensitive to the fine details of the momentum dis- 
tribution of the beam. We may take advantage of that 
by eliminating the specific dependence of the transverse 
wave numbers on the average beam radius I , and replace 
(ir /I) (n y , n z ) with the general symbol ~qtr, such that the 
specialized set of eigenfunctions, x n n z (yi z i x )i mav be 
replaced by a more general set X~q tT (y> z;x). 

The relativistic asymptotic (initial and final) energies 
of an electron 'trapped' within the beam double-cone 
boundary are: e? * = 



TOnC 4 + 



(fx? 



2 , h 2 TT 2 n 2 c 2 



2 4 , 

JflgC + 



W) 



W) +W) 



I 2 , , 
2 , where: 



W) 



(P 



The corresponding longitudinal momentum transfer is 
calculated from: 



[9x ~ 

(4 



4) 



(ql 



[qir 



j\ = 



where: 



(ei - £/) Oh + £/) / (tic) + 



(Qx + ql) ~ 2mv/h , (e, + £/) « (2mc 2 ), so that: 



\q x 
W, 

Mr 



ql) 2mv/h 



{e l -e s ){2mc 2 ) /{he) 2 + 



(q(r) - (qlY 



namely: 



Aq x w (lu/v) + h 



(q{r) - {OtrY 



/2mv 



VII. APPENDIX B 

In this appendix, following the method developed 
in ReL^ 2 -, we consider the dielectric loss function 



Im 



-2K*b 



K* 



appearing in Eq. (|lip . and show 

that it is proportional to KcE x (^K: —b\Luj- the electric 

field component along the e-beam axis at the beam posi- 
tion z = — b. The latter is the key ingredient in the cal- 
culation of the power loss function in the classical limit. 
We shall also show in this appendix that in the long 
wavelengths limit discussed in Ref.— the dielectric loss 
function reduces to the well known expression derived in 
Re&ii. 

Our analysis starts from the expectation value of the 
four- vector potential (A v = (cp, ~ A^J ), given by: 
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j,(odd,even) — ( 

2 



dd.even) 



A, t) = -^ [°° dt' f d 3 r'D^ (V , V; t - t') f xt ^ (V ', t') 
c u=o J— °° J 

^g-^ (odd,ev en) _ ^ -^j-(odd.euen) | j 

where j ext ^ ( r t') is the external four-current density 
generated by the e-beam (with the components j e:Et < 1/ = 

(c P '",7"<y, 



-1 



(B5) 



The definition of the matrix U can be found in RefJ^). 



I ), and D v>ll (r , r ';t — if) is the "dressed" In the limit of a semi-infinite slab (c* — > 00), r 

retarded photon Green's function, defined by the corre- r = r , so tnat. 
lator: 



(odd) _ 



(B2) 



(odd) , (etiem) r> * 



00 



The "bare" four- vector potential is given by: 
4 0) C?,i) = - [°° dt' [ d 3 r'D^ (I* -r^t-i') 



where D$ (V, V; t - f) = M° } (7* - V; i - i') £„ !jU , 
is the "bare" retarded photon propagator in the Lorentz 
gauge, which is given by: (V — ~~r'\ t — t') — 

rjj (t - if) 8{\~r -T'\/c + t- if) / [r - V\, with: 
1,^ = 



V* = 



-1 , v = 1,2,3 



The corresponding Fourier transforms with respect to 
the spatial coordinates parallel to the surface, with wave 
vector K — (k x ,k y ) are: 



A v [K- z;lu)=±J2J dz'D^ [K; z, z' 
Xj ext ^ (K;z';cj\ and 

,(0) 



D v 0) (K,z,z';u 



2-rK 



DV> (K*,\z-z'\) 

Our explicit expression for the external 4-current den- 
sity associated with the e-beam is: 



f xt >» (k;z';lo) 

-ce5 (z + b) 5 (vk x — to) 



(B3) 







z(f^5{z + b)8{vk x -uj), n = l 

, M = 2,3 



For z, z' < , i.e. both on the vacuum side of the 
dielectric slab, occupying the space: 2c* > z > , the 
lateral Fourier transform of Eq. (|B2p can be written in 
the form: 



K , z, z ; lo 



2irK* 



8v,u,e 



-K'\z-z' 



and: 



A v (K;z;u?j 



1 V„ 



z+b\ 



c2-kK* 
xceS (vk x — u>) 

( £ „-K*\z+h\ 



- 2r v ne 



K*(z-b) 



xe 



c2ttK 
lo 



6 (vk x — lj) 



- 27V ,1 e 



K*(z-b) 



where at z ■ 
A v (K--b-uu 



2ttK 



-5 (vk x - u) 7] u 



(B6) 



6, 
2rv n 



-2K*b 



(ri od n d) + rlT l) ) e K *( z+z ') J-ro,*% 
(B4) 



To simplify the calculation we shall consider in what 

follows the special case when the wave vector K is paral- 
lel to the c-bcam direction, which was selected along the 
x-axis, so that k y = , and K = k x . 

Thus, the electric field component along the e-beam 
axis at the beam position z = —b , E x (yK; —6; cuj , can be 

calculated from the explicit expressions for the potentials 
in Lorentz gauge, namely: 

E x (k- -b; w) = -ik x A (K; -b; lo] 

i(u/c)A! (K;-b;u?j 

-ik x (r ,o + (^) r ,i) 
1 (uj/c) (r 1)0 + (^) rx,! 

•0,0 + (^) r 0ll + (^) r li0 + (£) 2 n, 
or (by exploiting the symmetry property: ri.o 



-2K*b 



where the generalized reflection four- matrices for incident 
waves, which are either symmetric or antisymmetric with 
respect to the slab center, are given respectively by (see 
Rcf. 1 ^: 



E x ( K; —b; lo ) = — ik, 



ro,o 



LO 



r l,l 



~2K*b 



(B7) 
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This expression should be compared to the dielectric 
response function: 

r fA (k x ,k y ,u>) 

ro,o + (hql/mc) 7*0,1 + (hql/mc) r lfl 
+ [&<ll<fx/ ( mc ) 2 ) r hi 

which may be simplified (again due to the symmetry 
rifl = —r i and the inequality \Aq x \ as u>/v <C mc/h 
) to: 

1 x ' y ' ] " +{hAq x /mc)[r 0il -(v/c)r 1>1 ] 



rofl + (v/c) 2 ris 
+ (hu/mcv) [ro.i — (v/c) r% t 

ro,o + ( v /c) 2 r is 



(B8) 



Thus, since in the "classical" limit the prefactor of 7*1,1 
in Eq.(B7} < > {v/c) 2 , we find that: 



ReE x (K; —b;uj\ oc Im (7*0,0 + {v/c) 2 n,ij e 



2K*b 



Exploiting the continuity equation to connect various 
components of the matrix XJ^: 

k x U n - (uj/c)U i = -K*U 31 , and k x U 13 - 
(to/c) Uq3 = —K*U33, and noting the symmetry relation: 
U31 — —U13, it can be shown that: 



ro,o 



ro.i 



(uj/ck x ) U i + U { 



mi 



2(1 + TrU) 



Consequently: r 0fi + (v/c) 2 n,i -> r , + ^1,1 = 

2(l+TrU) 

and in the long wavelength limit, where 1 + 
TrU — ^ e ^ s x?^ and e is the bulk optical (frequency de- 
pendent) dielectric function of the platelet, we find: 



C0.0 



(1 



2s 



2K* 2 



(B9) 



^1, 



(s-1) (u/cY 



AK* 



1 + (e - 1) 



A" 



*2 



(0 + A*) z 



(BIO) 



so that finally: 



Re£ ; 



oc Im 



-2K'b 



K* 



K'(e-1) 



+ 



'v\ 2 ( (K*-Q) , (l~e)K 2 
yc) \(Q+K') (Q+K')(eK' + Q) 



which is equivalent to surface dielectric loss function 
obtained in RehW by using Maxwell's equations with 
macroscopic boundary conditions. 



Furthermore, the continuity equation also implies: 
k x U i + (uj/c) Uoo = K*U 30 , and k x U 13 - (u/c) U 03 
-K*Uz3, so that 
since C/30 = f/o3i we also find that: 



7*1,0 + 



U) 



Uqi - (m/ck x ) U 1 
2(1 + TrU) 
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